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Abstract 

Anderson localization of electron states on graphene lattice with diagonal and off-diagonal (OD) 
disorder in the absence of magnetic field is investigated by using the standard finite-size scaling 
analysis. In the presence of diagonal disorder all states are localized as predicted by the scaling 
theory for two-dimensional systems. In the case of OD disorder, the states at the Dirac point 
(E = 0) are shown to be delocalized due to the specific chiral symmetry, although other states 
(E ^ 0) are still localized. In OD disorder the conductance at ii^ = in an M x L rectangular 
system at the thermodynamical limit is calculated with the transfer-matrix technique for various 
values of ratio M/L and different types of distribution functions of the OD elements tnn'- It is 
found that if all the tnn'^s are positive the conductance is independent of L/M as restricted by 2 
delocalized channels at E = 0. If the distribution function includes the sign randomness of elements 
the conductivity, rather than the conductance, becomes L/M independent. The calculated 
value of the conductivity is around in consistence with the experiments. 
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I. INTRODUCTION 



As a potential candidate material for carbon-based electronics, graphene, a two- 
dimensional (2D) honeycomb lattice of carbon atoms, has attracted much attention 

nnnnn 

1|, 121, y, U, 15| . It is found that both the type (electrons or holes) and the number of carriers 
can be tuned by changing the gate voltage . Such tunability originates from the specific 
dispersion relation of the massless Dirac fermions in graphene |8|, [Oj. Its transport properties 
are very special. There exists a universal maximal resistivity for all samples with the Fermi 
evel near the Dirac point, independent of their shapes and mobility. The maximal resistivity 
sl implies the existence of channels of extended states at the thermodynamical limit, and 
the shape independence of the resistivity suggests that both the number and spatial config- 
uration of these extended channels depend on the shape. This behavior is certainly exotic 
as from the scaling theory all the states are localized in 2D systems [l^ . The non-zero mini- 
mal conductivity has been investigated by several theoretical studies using different methods 
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14 . Il5l |. The minimum conductivity was correctly demonstrated, but as to the 
reason for it and as to its exact value being 4e^//i or Ae'^/nh there is no consensus at present. 
Morozov et al. have found that in graphene the weak localization is strongly suppressed and 
they attributed this to mesoscopic corrugations of graphene sheets which can cause a de- 



phasing effect similar to that of a random magnetic field 



16j. Theoretically, the influence 



of disorder on 2D electron gases on the honeycomb lattice was studied with a self-consistent 



Born approximation JlTI, [l8| . Pereira et al. have investigated the localized states near va- 
cancies in graphene 19|. Khveshchenko has carried out an analysis of quantum interference 
effects in disordered graphene and specified the conditions for the quantum correction to 
the conductivity to be positive, negative, and zero [20I. By evaluating the dependence of 
the magnetoresistance of graphene on the relaxation rates from different sources McCann 
et al. showed that the warping disorder and the disorder due to atomically sharp scatterers 
have different effects on the magnetoresistance 21|. Morpurgo and Guinea investigated the 
effects of static disorder caused by curvature or topological defects and found that when the 
intervalley scattering time is long enough an effective time-reversal symmetry breaking can 
be induced 22| . Aleiner and Efetov showed two different types of logarithmic contributions 



to conductivity by describing the disorder with a renormalization-group equation 



231]. From 



a combination of mean-field and bosonization methods Altland pointed out that there are 
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several mechanisms conspiring to protect the conductivity although disordered graphene is 
subject to the common mechanisms of Anderson localization [OA]. From the linear response 
;heory Ziegler indicated that the transport properties in graphene are dominated by diffusion 



251]. Based on numerical study with finite-size Kubo formula for doped graphene, Nomura 



and MacDonald found distinct differences between the cases of short-range and Coulomb 
randomly distributed scatterers and speculated that the transport properties in graphene 
are related to the Boltzmann theory for Coulomb scatterers [26|. Sheng et al. carried out 
transfer matrix studies of graphene in a magnetic field to investigate the disorder effect and 
phase diagram [27|. Louis et al. calculated conductances for sarnples with edge disorder and 
obtained size independent conductivity as ~ Ae^/irh at E = 28|. Moreover, by adopting 
disorder model without intervalley scattering, Ostrovsky et al. constructed an effective field 
theory and showed that the system is at a quantum critical point with a universal value of 
the conductivity of the order of e'^/h [29]. With a similar model Titov demonstrated the 



3(3 ]. Bardarson et al. 



impurity-assisted tunneling by using the transfer matrix approach 
numerically showed that the conductivity at the Dirac point increases logarithmically with 
sample size in the absence of intervalley scattering [31]. From numerical calculation Nomura 
et al. demonstrated that all states are delocalized even in the strong disorder regime if the 
intervalley scattering is excluded [3^. San- Jose et al. investigated the full transport statis- 
tics of graphene at low dopings without intervalley scattering and found identical scaling of 
all current cumulants [33]. 

In this paper we investigate the Anderson localization of fermions at the Dirac point 
in graphene with different types of disorder by using the standard finite-scaling analysis 
together with the transfer-matrix technique. The Dirac point corresponds to the band cen- 
ter [E = 0) of a tight-binding Hamiltonian on the honeycomb lattice. In the real-space 
tight-binding representation the disorder can be introduced with randomly distributed site 
energies, as called diagonal disorder, or with random hopping integrals named as off-diagonal 
(OD) disorder. It is known that the localization behavior at the band center of a bipartite 
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361 ]. For instance. 



lattice with OD disorder is anomalous due to the chiral symmetry 
the decay of wavefunctions at the band center on 2D square lattice with OD disorder is 
much weaker than exponential. The honeycomb lattice is also bipartite, and there are other 
unusual features at the Dirac point such as the massless Dirac fermion dispersion relation. 
So it is interesting to investigate the difference in the localization behavior at the Dirac point 
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between diagonal and off-diagonal disorders. We derive the formula of transfer matrix for 
honeycomb lattice from which the numerical studies can be carried out to calculate inter- 
esting properties such as the rescaled localization length, the size and shape dependencies 
of the conductance. The finite-size scaling analysis reveals that although all the states are 
localized in the system with diagonal disorder, in the case of OD disorder, the states at the 
Dirac point are delocalized. From the calculated shape dependence of the conductance at 
= on an L X M rectangular lattice we find that when all the hopping integrals are 
positive, the conductance is independent of L/M, but if the distribution includes both the 
positive and negative hopping integrals, the conductivity becomes shape independent. Diag- 
onal disorder corresponds to atomic scatterers such as doped impurities, while OD disorder 
without sign changes of hopping integrals represents randomness of bond lengths caused 
by topological corrugations or warping in the sheets. For graphene OD disorder is more 
important than diagonal disorder owing to the ultrathin one-layer structure and the purity 
requirement in preparing the samples. In both types short-range disorder can exist and the 
intervalley scattering is not excluded. Thus the obtained results should be different from 
those by using models without intervalley scattering due to different symmetry classes. The 
sign randomness of hopping integrals can be produced by geometric or Berry phase ~ ±7r 
which is acquired by electron wavefunctions in every period of cyclic evolution around Dirac 
point in parameter space. The Dirac point can serve as the degenerate point for geometric 
phase in "slow" cyclic evolutions, such as in-plane vibrations of displacement between two 
sublattices, due to the specific dispersion relation and electron-lattice interaction. 



II. BASIC MODEL AND FORMALISM 



The main features of the band structure in graphene can be well described by a tight- 
binding Hamiltonian with one vr orbital per site on a honeycomb lattice sl: 

H = ^e„aTan ^ +tnn'{alan' + aj^/a„), (1) 

" {n,n') 

where aj^ (a„) creates (annihilates) an electron at site n, e„ is the energy level at site n, (. . .) 
denotes the nearest-neighbor (NN) sites, and tnn' is hopping integral between NN sites n 
and n'. Here, the spin indices are not explicitly included. In the absence of disorder, e„ = 
for all n and tnn' = to for all {nn') , this Hamiltonian leads to the 4-component Dirac fermion 



dispersion relation near the Dirac point E — 0. The diagonal and off-diagonal disorder can 
be introduced by adopting random variables e„ and tnn' satisfying distribution probabilities 
Prf(e„) and Poitnn'), respectively. These probability functions define strengths and types of 
disorder in the system. Diagonal disorder represents potential fluctuations due to impurities 
or due to randomness on substrate surface. A most commonly used distribution probability 
for diagonal disorder is the square function 



where W is the distribution width describing the strength of disorder. For OD disorder, there 
are different distribution functions describing different types of randomness. If we consider 
slight fluctuations of bond lengths around their average value due to lattice distortions, the 
following distribution is suitable 



where A (< 2|to|) is a measure of the bond- length fluctuations which make the NN hopping 
integrals randomly shift from to- 

In Poi{tnn') all hopping integrals have the same sign as to, reflecting that only changes 
of magnitude of tnn' a-re considered. In some cases, however, tnn' not only fluctuate in 
magnitude, but also change in sign due to speciflc physical mechanisms. One of these 
possible mechanisms is the geometric phase due to dynamical vibrations. The on-plane 
displacement between two sublattices, denoted as d = {dx,dy), changes all the NN bond 
lengths and thus linearly couples with electron states in the regime of small |d|. At a 
Dirac point {E — 0) the states are doubly degenerate at d = 0. The displacement d 7^ 
linearly removes the degeneracy, leading to conical energy surfaces in the parameter space 
spanned by and dy with the degenerate point situated at the origin d^ = dy = 0. Owing 
to thermal excitations or even to zero-point oscillations, ri^. and dy vibrate and both of 
them are dependent on time with period T = 2tt/uj where u is the frequency of vibrations. 
Usually the amplitudes and relative phase of dx and dy vibrations are statistically random. 
For generic values of them, the track in the dx-dy parameter space in a period is an ellipse 
encircling the degenerate point. Thus, according to the Berry theorem, in each period of 
vibrations electron wavefunctions at Dirac point will acquire a geometric phase ~ ±7r in 






= < 



1/A, for to - A/2 < tnn' <to + A/2, 
0, otherwise. 



(3) 
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addition to the ordinary dynamical phase. At the same time in this period the tunnehng 
electron has traversed a path of length vpT which equals vfT/Iq NN bonds with vp and /q 
being the Fermi velocity and bond length respectively. This means that the electron states 
change sign in traversing every vpT/Iq bonds. Such an effect occurs very naturally since the 
vibrations and electron transport are proceeding simultaneously. To include such an effect 
in a primary way, we introduce the sign-changing probability of an NN hopping s. This 
probability can be estimated from the length by which the electron state changes sign due 
to the geometric phase, i.e., s ~ Iq/{vfT). With such a sign-change model the wavefunction 
can automatically change sign in traversing every length of vpT which, averagely, contains 
only one hopping with opposite sign. This is a special gauge in which the geometric phase 
acquired in a path of vpT /Iq bonds is concentrated into one bond in this path. This gauge 
is valid under the condition that the interferences between paths shorter than vpT are not 
relevant. As we will see below, the transport properties at Dirac point in the case of OD 
disorder (including the sign changes) are scale-invariant. So the short-range processes are not 
important and the above condition can be satisfied. This sign-change mechanism is different 



from the dephasing effect in the model of Ref. [la where the distortions are fluctuating in 



the space (the ripples) but are static in time. We introduce the following distribution of tnn' 

Po{tnn') = sPol{tnn') + (1 " s)Pol ("^nn' ) , (4) 

to account for the sign-change probability s. When s = or s = 1, all the hopping integrals 
have the same sign and Poitnn') reduces to Poiiinn')- At s = 0.5 the probability of the sign 
changes is maximal. Below we set to as energy units. 
The wavefunction of a fermion can be written as 

^ = ^c„|n), (5) 

n 

where c„ is the amplitude on site n. At a given energy E the amplitudes are given by 
Schrodinger equation Hip = Eip. The propagation properties of an electron on a finite lattice 
can be given by a relation between amplitudes on both sides. This is usually formulated by 
the transfer-matrix method or by the Green function. We note that graphene lattice can 
be divided into two sublattices as distinguished by red and blue sites in Fig. 1. From the 
tight-binding nature of the Hamiltonian and the structure shown in Fig. 1 it can be seen 
that for given E all amplitudes can be iteratively calculated along the x {y) direction if the 
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FIG. 1: Structure of graphene lattice, i and j are indices of columns (shaded stripe along the y 
direction) and rows (shaded stripe along the x direction), and the overlapping area of two stripes is 
a basic cell where the lower and upper sites correspond to sublattices I = 1 and I = —1, respectively. 

amplitudes in the light-cyan shaded column (light-magenta shaded row) are given. From 
this a transfer matrix along the x or y direction can be given. We identify a site n with three 
integers, n = where i and j index the column and row, respectively, and / = ±1 

specifies the sublattice which the site n belongs to. In this notation the iterative relation 
for amplitudes along the x direction can be derived from Schrodinger equation as 



For a system with finite width in the y direction, one can define vectors for columns by 
Uj ^ = (ci,i,z, Q,2,h • • • , Ci^M,iV ) where M is the number of cells in a column. Then Eqs. ([6]) 




(6) 



Ci+l,j,l 



{E — ej_(_ij^_i)cj+ij^_i — tj+ij,-i;ij,iCij^i — tj+ij^_i;j_|_ij+(_i)i^iCj_(_ij+(_i)i_i 



(7) 
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and ([7j) can be written in a matrix form: 

Ui,i = Ta^jUj+i^i + T4^jUj_|_i^_i, 
where elements of matrices are given by 



t 



ti,j,l;i+l,j,-l 
ti+l,j,-l;i+l,j,l 



k,j,l;i,j + {-iy 



ti,j,l;i+l,j,-l 



-6 



i+ l,j-l;i+l,j + {-iy,l ^ 



(8) 
(9) 

(10) 

(11) 
(12) 

(13) 



Here, one may use the periodic or open boundary condition at the ends j = 1 and j = M. 
We have the iterative relation for whole column 



(14) 



1 1 













where the transfer matrix is 

t = 



V 



(15) 



T} — \r-p rp rri — 1 rji — lrji rji 

If there are L columns in the x direction, the total transfer matrix which gives the relation 
between amplitudes at two end columns is 



f =nTL-.+i- 

2=1 



(16) 



If we want to investigate the propagation properties in the y direction, we have to establish 
the transfer matrix for rows. The relation between two adjacent rows is written as 



(17) 



where Vj_;'s are vectors denoting rows with length M, Vj^^ = (cij,;, C2jy, . . . , cmj,/)"^, and the 
elements of the transfer matrix are 



;i8) 



where 



t 



a 



i,j,{-iy;i,j+i-(-iy 



(2) 



t 



c 



(4) 



E 



t 



-(-1)" 
k,j+i,-{-iy;i,j,{-^y 



^(3) 



t 



i,j,{''^y;i+{~iy,j-i-iy 



ti 



c 



(6) 



The total transfer matrix for L rows is 



V 



i=l 



(19) 



For a finite M x L system, there are M propagating channels along the length L. The 



Hermitian matrix, (Ttj')i/2i qj. ^yty^i/2L^ ^i^s 2M eigenvalues. These 2M eigenva. 



ues 



should be positive and come in inverse pairs due to the unitarity of the matrix [37|, |38|, |39|, |40| . 
The M positive logarithms of eigenvalues, denoted as 7^ in an ascending order with index i, 
are Lyapunov exponents (LE) reflecting the exponential decay of the corresponding channels. 
Obviously, the logarithms of other M eigenvalues are the negatives of 7j. The LEs will 
approach when the disorder vanishes. The localization length is defined as the reciprocal 
of the smallest LE ^ = I/71. The localization length plays a crucial role in the localization 
theory. 

The zero-temperature conductance (7 of a rectangular MxL system in units of 2e^/ h, with 
;wo spins being taken into account, can be evaluated using M-channel Landauer formula 
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43| 



9 = MiH), (20) 

where the M x M transmission matrix t describes the transmission of electrons from one 
lead to the other in the longitudinal direction. Here we assume that the leads are connected 
to all the M channels and they are perfect metals with band width much larger than that 

nn 

of graphene. Formula (120]) can be expressed with the LEs as pTj, 4Q[ 



M 



(21) 



cosh^(7iL) 

Numerically, LEs can be calculated by using the standard method of Gram-Schmidt 
re-orthonormalization after each a few, say, ten steps of multiplication of transfer matrices 
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45 1 ■ This is equivalent to the diagonahzation of T^T or V'^V, but avoids terrible overflow 



and loss of precision on the computer. On the other hand, although the calculations of 
transfer matrices along the x and y directions are certainly different, there is no essential 
difference in the calculated results between these two directions. This is owing to the fact 
that the dispersion relation is conic around the Dirac point which is isotropic in the x and 
y directions. Below we will only present the results obtained for the transmission along the 
y direction. 



III. LOCALIZATION IN SYSTEM WITH DIAGONAL DISORDER 

First we investigate the case of diagonal disorder, i.e., W ^ 0, X = 0, and s = 0. The 
lowest LE 7i(M) is calculated for very long (L ~ 10^) strip of width M. The rescaled 
localization length is deflned as ^^^^ = For given energy and disorder strength, 

the localization behavior of the fermions can be determined by the scaling dependence 
of If ^^^^ increases with M, the states are extended. For localized states ^^^^^ is 

decreased by increasing M. The M-independence of ^^^^ corresponds to the flxed point 
separating regions of extended and localized states. In Fig. 2(a), we plot the M dependence 
of ^^^p for states at = in systems with different strengths of diagonal disorder. For all 
investigated strengths of diagonal disorder, the rescaled localization length always decreases 
with increasing width. The large fluctuations in the curve of W = 1 are due to numerical 
errors in the Gram-Schmidt re-orthonormalization procedure which become serious in the 
case with ^ > 1. This means that the states at the Dirac point are localized by tiny 
diagonal disorder. In Fig. 2(b), we plot the M-dependence of ^^^^ for a given strength 
of diagonal disorder but different energies. Although all states are localized, the scaling 
behavior is much different between the states at Dirac point and the other states: the 
decreasing slope of the former is much larger than that of the latter. This implies that the 
states at Dirac point are more easily localized than other states by introducing the diagonal 
disorder in the thermodynamical limit, in spite of their larger localization length in smaller 
systems. This may be due to the small density of states at E = and the breaking of the 
chiral symmetry by diagonal disorder. The localization behavior of the states is consistent 
with the common theories of the Anderson localization in 2D. 
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FIG. 2: (Color online) Scaling behavior of rescaled localization length in the presence of diagonal 
disorder, (a) versus M at the Dirac point (E = 0) for different strengths of disorder, (b) 

^^"^^ versus M for a given strength of disorder and various energies. 

IV. DELOCALIZATION AT DIRAC POINT IN THE CASE OF OFF-DIAGONAL 
DISORDER WITHOUT SIGN CHANGES 

Differently from diagonal disorder, off-diagonal disorder does not break the chiral symme- 
try at Dirac point, leading to an essential difference in the localization. The scaling analysis 
of the rescaled localization length and the conductance at Dirac point is shown in Fig. 3 
by using the distribution function Poi in Eq. ([3]). This corresponds to the case of randomly 
warped or corrugated lattice where bond lengths are randomly shifted from the perfect value. 
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From Fig. 3(a) we can see that the rescaled locahzation length is independent of M for all 
disorder strengths. This implies that the localization length at the thermodynamical limit 
is infinite, corresponding to delocalized states at = 0. The conductance of an M x M 
system calculated from Eq. fl?T]) is also independent of M for all investigated strengths of 
off-diagonal disorder [see Fig. 3(b)]. So the conductance at the thermodynamical limit can 
be determined from these size independent values. The dependence of g on the disorder 
strength is shown in the inset of Fig. 3(b). For relatively weak off-diagonal disorder (A < 1) 
the conductance at Dirac point is about corresponding to the quantum conductance of 
4 channels (2 spin channels and 2 equivalent channels from the chiral symmetry at = 0). 
In fact, in a system with off-diagonal disorder, the calculated Lyapunov exponents at = 
can be grouped into pairs, i.e., 72i-i = 72i for integer i in the range between 1 and cor- 
responding to two equivalent channels with the chiral transformation. These channels have 
the localization length proportional to the width M. At the thermodynamical limit they 
become infinite. By increasing the disorder strength from A = 1 towards its maximum value 
A = 2, these channels are still delocalized at the thermodynamical limit but the conductance 
rapidly drops, as can be seen from the inset of Fig. 3(b). 

The delocalization behavior at = in graphene with the off-diagonal disorder is not 
consistent with the usual localization behavior of 2D disordered systems. One may attribute 
this inconsistency to the chiral symmetry at E = Q which is conserved in off-diagonal disorder 
but destroyed in diagonal disorder. So the localization at energies 7^ 0, where the chiral 
symmetry no longer exists, could be expected. In Fig. 4, we plot the scaling behavior 
of the rescaled localization length and the conductance of an M x M system for different 
energies. Except for the curves of = 0, both the rescaled localization length and the 
conductance decrease by increasing M in the general trend if fluctuations of some curves are 
ignored, suggesting the usual localization behavior aX E ^ From Fig. 4(a), the rescaled 
localization length is generally decreased by increasing E from = for a fixed M. This 
means that in off-diagonal disorder the states of energies farther apart from Dirac point 
are more localized. From Fig. 4(b), however, for a given value of M the conductance at 
larger \E\ has larger values. This can be explained by the fact that there are more electron 
channels at larger ji^'l in graphene which can contribute to the conductance if M is not very 
large. 

By using a scaling transformation M —>■ M/^ with appropriate value of scaling factor ^ 
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FIG. 3: (Color online) Scaling behavior of the rescaled localization length (a) and the conductance 
of an M X M system (b) at the Dirac point {E = 0) in the presence of off-diagonal disorder. The 
inset of (b) shows the dependence of g on the strength of off-diagonal disorder A. 

for every energy E, the curves in Fig. 4(a) for different energies can merge together. In Fig. 
5 we plot the common dependence of ^^^^ on M/^ for various energies. The scaling factor 
^ as a function of E is shown in the inset which can be fitted by a simple function ,^ = ^ 
with prefactor a = 2.66. Since ^ is just the localization length at the thermodynamical 
limit, this scaling behavior further ensures the delocalization a.t E = and describes how 
the localization length diverges when E 0. From the Dirac dispersion relation the wave 
length of electrons is also inversely proportional to \E\. This suggests that in the case of 
off-diagonal disorder the localization length and the electron wave length have the same E 
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FIG. 4: (Color online) Scaling behavior of the rescaled localization length (a) and the conductance 
of an M X M system (b) for different energies in the presence of off-diagonal disorder. 



dependence and can be regarded as the same length scale. We also notice that although the 



states at E = are delocalized, there is no real metal 
is certainly different from the results in Refs. 
free of intervalley scattering are adopted. 
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ic region with finite range of E. This 
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where the disorder models 



The different localization behaviors between E = and E ^ can be more clearly 
depicted by the spatial distributions of wavefunctions as shown in Fig. 6. As can be seen 
from Fig. 6(a), the delocalized state at = is more likely to be a Cantor-set like critical 
state, rather than a real metallic state. Such states may provide propagation paths for 
current from one edge to the other. On the other hand, the state at E = 1.8 shown in 
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FIG. 5: (Color online) Common dependence of the rescaled localization length on the scaled 
width f for various energies. Inset: ^ as a function of E (circle symbols) which is fitted by function 
^ = (red curve). 

Fig. 6(b) has significant amplitudes only in some isolated islands, corresponding to typical 
localization. 

In the case of off-diagonal disorder there are 4 delocalized channels at Dirac point, giving 
conductance This leads to almost shape-independent conductance and shape-dependent 

attice investigated by Tworzydlo 



conductivity. The situation is similar to the case of perfect 



et al. in Ref. 
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we calculate the conductivity 



46|. Following the procedure proposed in Ref. 
from the conductance of a M x L rectangular sample. The obtained shape dependence of the 
conductivity is shown in Fig. 7. This behavior is consistent with the nature of the quantum 
transport in 2D, but is inconsistent with the experiments where the conductivity is shape 
independent. This means that to explain this experimental result we have to invoke some 
mechanism which could convert the quantum transport to the classical-like transport. 



V. CLASSICAL-LIKE TRANSPORT AT DIRAC POINT IN THE CASE OF OFF- 
DIAGONAL DISORDER WITH SIGN CHANGES 

The quantum transport processes could be drastically changed if the phases of wavefunc- 
tions become uncertain. This may happen due to the geometric phase which can be acquired 
by electrons during a cyclic evolution of "slow" degrees of freedom around a conic point. 
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FIG. 6: (Color online) Contour plot of spatial distributions of wavefunctions at (a) E = and (b) 
E = 1.8 in a 60 X 60 system in off-diagonal disorder of strength A = 0.7. 

The additional geometric phase of ~ ±7r corresponds to sign changes of the wavefunction 
during the motion of tunneUng electrons. In the diffusion regime, if the sign changes ran- 
domly happen during the spatial motion, they can be modeled by a random distribution 
of positive and negative signs of hopping integrals. In probability Po{tnn') of Eq. (jlj) we 
include both the sign randomness and the magnitude randomness of the hopping integrals 
which are specified by parameters s and A, respectively. 

First we set A = and focus on the effect of sign randomness. In Fig. 8 we plot the 
scaling behavior of the rescaled localization length and the conductance of an M x M system 
for different values of s. Similarly to the case of off-diagonal disorder, the states at the Dirac 
point in the presence of the sign randomness are also delocalized in the sense of finite-size 
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FIG. 7: (Color online) Width dependence of conductivity of an M x L rectangular system with a 
given length at E = for different strengths of off-diagonal disorder. 

scaling analysis. Differently from the behavior in the off-diagonal disorder shown in Fig. 3, 
effects of sign randomness on the rescaled localization length and on the conductance are 
opposite: ^^^P in s = 0.5 is smaller than that in s = 0.2, but g in s = 0.5 is larger than 



that in s = 0.2. This is strange because ^^^^ corresponds to the channel with the largest 
localization length which usually provides the leading contribution to the conductance. The 
only explanation for this is that the other channels also give significant contributions in the 
case of sign randomness. An expectable consequence of this is a new shape dependence of 
the conductivity, since the number of the tunneling channels is no longer restricted. 

The shape dependence of the conductivity can be investigated by calculating the conduc- 
tance of M X L rectangular samples. The obtained conductivity as a function of ratio M/ L 
for various sizes is shown in Fig. 9. Except for small values of M/L, the conductivity has 
values near independent of M/L and M, in consistence with the experirneiital findings. 



Surprisingly, this value is even larger than ^ obtained in ballistic graphene |46|, |47| in spite 
of the disorder introduced by the sign randomness of hopping integrals. To get insight into 
this anomalous feature, we adopt the average spacing between successive pairs of Lyapunov 
exponents denoted as ^7, i.e., 721+1 — 72i-i ~ 72i+2 — l2i ~ ^1- Notice that in this case the 
Lyapunov exponents are grouped into pairs, corresponding to two equivalent channels from 
the chiral symmetry. By keeping only the leading exponential term of cosh^(7jL), Eq. (l2Tl) 
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FIG. 8: (Color online) Scaling behavior of the rescaled localization length and the conductance 
of an M X M system in the presence of sign randomness of hopping integrals obeying probability 
Po{tnn') with different values of s. 

can be approximately rewritten as 

g « 8e--" « " K (22) 

1=1 ^ 



in units of In a perfect system of width M and using the periodic boundary condition, 
the spacing ^7 for states at = is ^7 = 27r/M due to the Dirac fermion dispersion relation 
One may reasonably suppose that ^7 oc 1/M is still held in the case of OD disorder. 
Then, for M ^ L, the conductivity of an M x L rectangular sample can be calculated as 



a 



M 6'yM 



(23) 



So for the perfect lattice o" = - in units of ^ |46|,|47|. In the case of sign randomness, the 
values of MSj are in the range [0.9, 2.5] as shown in the inset of Fig. 9. Together with the 
factor e~'^'^^^, the obtained conductivity may be vr times of that of the perfect lattice. Thus, 
the 7T factor difference in the conductivity is not trivial and reflects an essential change of 
the fermion properties due to the phase uncertainty. 

It is interesting to investigate how the localization behavior changes in the presence 
of the magnitude randomness of hopping integrals. The scaling behaviors of the rescaled 
localization length and the conductance of an M x M system with both nonzero s and A 
are shown in Fig. 10. It can be seen that the states at ii^ = are still delocalized and 
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FIG. 9: (Color online) Calculated conductivity as a function of ratio M jL for rectangular sheets 
with different widths M. Inset in the lower panel: Distribution function P{M6^) of M5^. The 
green dashed line indicates the 5-function distribution at 27r in the case of perfect sheet. 

the opposite effects on the rescaled locahzation length and on the conductance are more 
evident than those shown in Fig. 8. In Fig. 11 we display the shape dependence of the 
conductivity of the M x L rectangular samples for s = 0.5 and different nonzero values of A. 
Similarly to the case of A = 0, the conductivity is still shape independent except for small 
values of M/L. This means that the sign-randomness induced shape independence of the 
conductivity is robust against the warping or corrugation disorder which causes randomness 
of magnitudes of hopping integrals. It can also be seen that the value of the conductivity is 
slightly reduced from ^ by increasing the warping or corrugation disorder. For reasonable 
strength of the warping disorder (A < 1), however, the conductivity is above in the 
range of measured values in experiments. 

It is also interesting to investigate the scaling behavior of the statistics of the minimum 
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FIG. 10: (Color online) Scaling behavior of the rescaled localization length and the conductance of 
an M X M system in the presence of both sign randomness and magnitude randomness of hopping 
integrals obeying probability Po(t„„/). 




M/L 



FIG. 11: (Color online) Calculated conductivity as a function of ratio W/L for rectangular sheets 
with different widths M in the presence of both sign randomness and magnitude randomness of 
hopping integrals. 

conductivity. In Fig. 12 we plot the statistical distribution of the conductivity at £" = for 
various system sizes and for given values of M/L, A, and s. We note that the distribution 
of the conductivity is Gaussian-like, and both the average and the distribution width are 
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FIG. 12: (Color online) Statistical distribution of the conductivity at E = for various system 
sizes in OD disorder with sign randomness. The values of M/L, A, and s are given in the figure. 
The statistics is taken from 500 random configurations. 

almost scaling invariant. This is different from the statistics of a single Lyapunov exponent in 
usual disordered system for which the distribution width is 1/M dependent. This difference 
originates from the fact that in the present system the conductivity is contributed not only 
from the channel with the smallest LE but also from other channels. As a consequence the 
statistics of the conductivity given by Eq. fl23l) is mainly determined by the distribution 
of 6'jM rather than by the distribution of a single LE. It can be easily verified that in the 
present case the series of successional LEs times M, 71M, 72M, . . . , are scaling invariant. 
Thus the statistics of 57M is also scaling invariant, leading to the scaling invariant statistics 
of the minimum conductivity. 



VI. CONCLUSIONS 



We investigate the effect of different types of disorder on the transport properties of elec- 
trons at Dirac point in graphene. From the finite-size scaling analysis together with the 
transfer-matrix technique, it is found that all states are localized in the presence of the diag- 
onal disorder, in consistence with the main conclusions of the scaling theory for 2D. In the 
case of off-diagonal disorder, however, the states at ii^ = are delocalized. This delocaliza- 
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tion, inconsistent with the scahng theory and only existing at Dirac point {E — 0), is mainly 
due to the specific chiral symmetry which is conserved in off-diagonal disorder but destroyed 
by diagonal disorder. The obtained conductance at the thermodynamical limit shows two 
transparent chirally symmetric quantum channels at E = 0. These channels, however, could 
not produce shape independent conductivity observed in experiments. It can be obtained 
only by introducing the sign randomness of hopping integrals. We investigate the scaling 
behavior of the rescaled locahzation length and the conductance in the presence of both 
sign randomness and magnitude randomness to show the existence of delocalized states at 
E — 0. Especially, the calculated conductivity of a rectangular sample is shape independent. 
As the sign randomness of the hopping integrals corresponds to some additional uncertain 
phases, the obtained results reflect a general trend of transition from quantum to classical 
transport in dephasing processes. Although the calculated conductivity is in consistence 
with the experimental findings, the mechanism for the phase uncertainty still needs to be 
identified by future investigations. 
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